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We present a measurement of the Cabibbo-Kobayashi-Maskawa matrix element \V c b\ and the 
form-factor slope p 2 in B — > D£~ ' vt decays based on 460 million BB events recorded at the T(4S) 
resonance with the BABAR detector. B — > DC vt decays are selected in events in which a hadronic 
decay of the second B meson is fully reconstructed. We measure the differential decay rate and 
determine £(l)|V c6 | = (43.0 ± 1.9 ± 1.4) x 10~ 3 and p 2 = 1.20 ± 0.09 ± 0.04, where G(l) is the 
hadronic form factor at the point of zero recoil. We also determine the exclusive branching fractions 
and find B(B~ -> D a Tvi) = (2.31 ± 0.08 ± 0.09)% and B(B° -> D + i~V t ) = (2.23 ± 0.11 ± 0.11)%. 
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In the Standard Model (SM) of electroweak interac- 
tions, the rate of the semileptonic B — > D£~vt decay 
is proportional to the square of the Cabibbo-Kobayashi- 
Maskawa (CKM) Q matrix element | V c b | , which is a mea- 
sure of the weak interaction coupling of the b to the c 
quark. The length of the side of the unitary triangle op- 
posite to the well-measured angle (3 is proportional to the 
ratio |Kif,/V^b|, making the determination of |V c t>| impor- 
tant to test the SM description of CP symmetry violation. 
In addition, imprecise knowledge of \V c b\ is an important 
uncertainty limiting comparison of measurements of CP 
violation in K meson decays with those in B meson de- 
cays @| . 

Measurements of \ V C J have been performed in inclusive 
semileptonic B decays Q and in the exclusive transitions 
B — > D£~v t and B — * D*£~v g The most recent in- 
clusive and exclusive determinations differ by more than 
two standard deviations, with the inclusive result more 
than twice as precise as the exclusive one Q. Thus im- 
provements in the measurements of exclusive decay rates 
are highly desirable, particularly for B — ► D£~vg decays, 
where the experimental uncertainties dominate. Mea- 
surements of \V c b\ based on studies of B — > D£~Dg decays 
have previously been reported by the Belle [|| , CLEO Q , 
ALEPH H and BABAR 4] Collaborations. 

The B -> Dl'vt decay rate @ is proportional to the 
square of the hadronic matrix element that describes the 
effects of strong interactions in B — > D transitions. In 
the limit of very small lepton masses (£ = e or fi), their 
effect can be parameterized by a single form factor Q (w) : 

dw 

= ^ M UMb + M D ) 2 (w 2 - if' 2 I V cb | 2 G 2 (w), 



where Gf is the Fermi coupling constant, and Mb and 
Mb are the masses of the B and D mesons, respec- 
tively. The variable w denotes the product of the B and 
D meson four- velocities Vb and Vd, w — Vb • Vd — 
{M% + M 2 D - q 2 )/(2M B M D ), where q 2 = (p B - p D ) 2 , 
and pb and po are the four- momenta of the B and D 
mesons. 

In the limit of infinite quark masses, G(w) coincides 
with the Isgur-Wise function This function is nor- 
malized to unity at zero recoil, where q 2 is maximum. 
Corrections to the heavy quark limit have been calculated 
based on unquenched and quenched lattice QCD [l^] • 
Thus I V c b\ can be extracted by extrapolating the differen- 
tial decay rate to w = 1. To reduce the uncertainties as- 
sociated with this extrapolation, constraints on the shape 
of the form factor are necessary. Several functional forms 
have been proposed 1 1 31) . We adopt the parameterization 
suggested in Ref. [14J, where the non- linear dependence 



of the form factor on w is expressed in terms of a single 
shape parameter, the form- factor slope p 2 . 

In this letter, we present a measurement of the differ- 
ential decay rates for B° — + D + £~Dg and B~ — > D°£~P{, 
decays. The analysis is based on data collected with the 
BABAR detector [15J at the PEP-II asymmetric-energy 
e + e~ storage rings. The data consist of 417 fb _1 recorded 
at the T(4S') resonance, corresponding to approximately 
460 million BB pairs. An additional sample of 40 fb _1 , 
collected at a ccntcr-of-mass (CM) energy 40 MeV below 
the T(AS) resonance, is used to study background from 
e + e~ —*//(/ = u, d, s, c, t) continuum events. We also 
use samples of GEANT4 Monte Carlo (MC) simulated 
events that correspond to about three times the data 
sample size. The simulation models B — > D^i^De de- 
cays using calculations based on Heavy Quark Effective 
Theory (HQET) Q, -B ^ D**{^ D^-k)1~vi decays 
using the ISGW2 model and non-resonant B — » 

D^tt£~D£ decays using the Goity-Roberts model [l7j |. 
The MC simulation includes radiative effects such as 
bremsstrahlung in the detector material. QED final-state 
radiation is modeled by PHOTOS []J]. 

Semileptonic decays are selected in BB events in which 
a hadronic decay of the second B meson (B tag ) is fully 
reconstructed. This leads to a very clean sample of 
events and also provides a precise measurement of q 2 
and hence w. We first reconstruct the semileptonic B de- 
cay, selecting a lepton with momentum in the CM frame 
p\ larger than 0.6 GeV. For electrons, we search for a 
vertex formed in conjunction with a track of opposite 
charge and remove those with an invariant mass con- 
sistent with a photon conversion or a ir° Dalitz decay. 
Candidate D° mesons that have the correct charge cor- 
relation with the lepton are reconstructed in the K~ir + , 
K-iT+ir , K-ir+ir+ir-, ^Tr+Tr", K^ir+n-ir , K%ir°, 
K + K~ , 7r + 7r~, and KgKg channels, and D + mesons in 

the K-TT+TT+, K-TT+1T+1T , K%1T + , J^TT+TT , K+R-TT+, 

K° S K+, and K^tt+tt+tt- channels. D°(D+) candidates 
are required to be within 2a of the nominal D°(D + ) 
mass, where a is approximately 8 MeV/c 2 . In events 
with multiple D£~ combinations, the candidate with the 
largest D£~ vertex probability is selected. Events with 
more than one reconstructed lepton with p\ > 0.6 GeV 
are vetoed. 

We reconstruct B tag decays [l|| in charmed hadronic 
modes B — > D^Y, where Y represents a collection of 
hadrons, composed of niir + 71471" , 

where n\ + n 2 — 1,3,5, 773 < 2, and 774 < 2. Using 
D a (D+) and D*°(D*+) as seeds for B~(B ) decays, we 
reconstruct about 1000 different decay modes. The kine- 
matic consistency of a -B tag candidate with a B meson 
decay is evaluated using two variables: the beam-energy 
substituted mass m,ES = ^/s/4 — (pjj) 2 , and the energy 
difference AE = E* B — y/a/2. Here y/s is the total CM en- 
ergy, and p* B and E* B denote the momentum and energy 
of the Btag candidate in the CM frame. For correctly 



5 



identified B tag decays, ities peaks at the B-meson mass, 
while AE is consistent with zero. We select B t8ug candi- 
dates in the signal region defined as 5.27 GeV < m ES < 
5.29 GeV, excluding those that have daughters in com- 
mon with the B — > D£~i>£ decay. In the case of multiple 
Stag candidates in an event, we select the one with the 
smallest \AE\ value. The B tag and the D£~ candidates 
are required to have the correct charge-flavor correlation. 
Cross-feed events, i.e., S^ g (S° ag ) candidates erroneously 
reconstructed as a neutral (charged) B meson, are sub- 
tracted using estimates from the simulation. 

Semileptonic B decays are identified by 
their missing-mass squared value, m miss = 
\p(T(AS)) -p(B Us ) -p{D) -p(l)] 2 , defined in terms 
of the measured particle four-momenta. For correctly 
reconstructed signal events, the only missing particle is 
the neutrino and w^iss peaks at zero. Other semilep- 
tonic B decays, like B — > D^*'**} £~ where at least one 
particle is not reconstructed (feed-down), yield larger 
values of m^ iss . 

We measure \V c b\ and the form- factor slope p 2 by a 
fit to the w distribution. We examine the data and MC 
events in ten equal-size w bins in the interval 1 < w < 
1.6. Since the B momentum is known from the fully 
reconstructed B tag in the same event, w can be recon- 
structed with good precision, namely to ~ 0.01, which 
corresponds to about 2% of the full kinematic range. 

To obtain the B — > D£~V£ signal yield in each bin of 
w, we perform a one-dimensional extended binned maxi- 
mum likelihood fit [13] to the corresponding m 2 niss distri- 
bution. The fitted data samples are assumed to contain 
four different types of events: B — > D£~vi signal events, 
feed-down from other semileptonic B decays, combina- 
torial BB and continuum background, and fake lepton 
events (mostly from hadronic B decays with hadrons 
misidentified as leptons). The Probability Density Func- 
tions (PDFs) are derived from the MC predictions for 
the different semileptonic B decay m^ iss distributions. 
We use the off-peak data to provide the continuum back- 
ground normalization. The shape of the continuum back- 
ground distribution predicted by the MC simulation is 
consistent with that obtained from the off-peak data. 

The measured m^ iss distributions are compared with 
the results of the fits for two different w intervals in Fig. 1. 

We perform a least-squares fit to the observed signal 
yields in the ten bins of w. We minimize a \ 2 defined as 

X 2 = E " V ' -y ( 2 ) 

where the index i refers to the w bin and the index j 
runs over all MC events in bin i; iV| ata is the observed 
number of signal events found in the i th bin and <r^ ata the 
corresponding uncertainty. The expected signal yields 
are calculated at each step of the minimization from the 
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FIG. 1: (Color online) Fit to the m^ iss distribution, in two 
different w intervals, for B~ — > D°l~vf. the data (points 
with error bars) are compared to the results of the overall fit 
(sum of the solid histograms). The PDFs for the different fit 
components are stacked in the order shown in the legend. 

reweighted sum of jVj^ c simulated events. Each weight 
is the product of two terms, Wj — W c x Wj' thco where 
W is an overall fixed scale factor, which accounts for 
the relative integrated luminosity of the data and signal 
MC events, and Wj' 1 eo is computed using the true w 
value of the event j and depends on £7(l)|Vd>| and p 2 , 
which are free parameters determined in the fit that are 
recalculated at each step of the minimization. 

We first fit the w distributions for the charged and 
neutral B — > D£~vi samples separately, and then per- 
form a fit to the combined B — > D£~vg sample. In Fig. 
[2] we show the comparison between the data and the fit 
results for the combined sample. The measured values 
of Q{1) | V c b | and p 2 , with the corresponding correlation 
Peon obtained from the fit, are reported in Table |TJ The 
value of the branching fraction is computed by integrat- 
ing the differential expression in Eq. [1] and dividing by 
the appropriate _B-meson lifetime. 

In order to reduce the systematic uncertainty on the 
measurement of £?(1)|V^>| and the branching fractions, we 
normalize the exclusive signal yield to the yield of inclu- 
sive semileptonic decays, B — > X£~Di, in events tagged 
by a fully reconstructed hadronic B decay. The inclu- 
sive B — > Xl~vi decays are selected by identifying one 
charged lepton with p\ > 0.6 GeV and the charge ex- 
pected based on the B ta . g decay. In the case of multiple 
Stag candidates in an event, we select the decay mode 
with the highest purity, estimated from the MC predic- 
tion for the fraction of true decays in the m es signal re- 
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TABLE I: Fit results for each sample. In the last column we report the results for the B° and B~ combined fit, where the 
branching fraction refers to B° decays. We also report the signal yields and the reconstruction efficiencies, integrated over the 
full w range. For (7(l)|V" c t>| and p 2 , we report both the statistical and systematic uncertainties. 





B~ -» D a £-v t 


B° D+Tvi 


B -» D£-&£ 


0(1)|V,*| • 10 d 


41.7±2.1 ±1.3 


45.6± 3.3±1.6 


43.0± 1.9±1.4 


P 2 


1.14± 0.11±0.04 


1.29± 0.14±0.05 


1.20± 0.09±0.04 


Pcorr 


0.943 


0.950 


0.952 


X 2 /ndf 


3.4/8 


5.6/8 


9.9/18 


Signal Yield 
Recon. efficiency 
B 


2147 ± 69 1108 ± 45 
(1.99 ± 0.02) x 10~ 4 (1.09 ± 0.02) x 10~ 4 
(2.31± 0.08 ± 0.09)% (2.23± 0.11 ± 0.11)% 


(2.17± 0.06 ± 0.09)% 




FIG. 2: (a) w distribution obtained summing B~ — > D°£~i>i 
and B° — > D + £~vi yields. The data (points with error bars) 
are compared to the results of the overall fit (solid histogram), 
(b) £?(w)|V c i,| distribution corrected for the reconstruction ef- 
ficiency, with the fit result superimposed. 



gion. Background components that peak in the tues sig- 
nal region include cascade B meson decays, for which the 
lcpton does not come directly from the B, and hadronic 
decays, and are subtracted using the corresponding MC 
predictions. The inclusive B — » X£~i>t yield is obtained 
from a maximum likelihood fit to the tues distribution 
of the _B tag candidates, as described in Ref. [5l|]. The fit 
yields (198.9±1.6)x 10 3 events for the B~ -> X£^v t sam- 
ple and (116.3 ± 1.0) x 10 3 events for the B° -> X£-y t 
sample. The corresponding reconstruction efficiencies, 
including the B tag reconstruction, are 0.39% and 0.25%, 
respectively. 

Numerous sources of systematic uncertainty are inves- 
tigated. The largest uncertainties are due to differences 
in the efficiency of the B tag selection between the ex- 
clusive B — > D£~V£ and inclusive B — > X£~vi decays 



(a relative 1.5% systematic uncertainty on \V c b\), the 
B — > D£~i>i fit procedure (1.3%), and the uncertainties 
on the branching fractions of the reconstructed D decay 
modes and B — > D**£~D e decays (1.1%). The uncer- 
tainties due to the detector simulation are established 
by varying, within bounds given by data control sam- 
ples, the tracking efficiency of all charged tracks (0.7%), 
the calorimeter efficiency (0.9%), and the lepton iden- 
tification efficiency (0.9%). We evaluate the systematic 
uncertainties associated with the MC simulation of var- 
ious signal and background processes: photon conver- 
sion and 7T° Dalitz decay, B cascade decay contamination 
(0.8%), and flavor cross-feed (0.2%). The uncertainty 
arising from radiative corrections (0.1%) is studied by 
comparing the standard results with those obtained when 
PHOTOS is not used. We take 30% of the difference 
as a conservative systematic uncertainty. We vary the 
B — > D*£~T>i form factors (0.4%) within their measured 
uncertainties Q and use an HQET parameterization 22 1 
to describe B —> D**£~vn decays (0.3%). We evaluate 
an uncertainty associated with the B — > Xl~vi fitting 
procedure (0.8%), and with the absolute branching frac- 
tion B(B — * Xl~v£) used for the normalization (0.8%). 
The complete set of systematic uncertainties is given in 
Ref. [23. 

From the fit to the combined B — > D£~vi sample, 
we measure Q(\)\V ch \ = (43.0 ± 1.9 ± 1.4) x 10" 3 . Us- 
ing an unquenched lattice calculation (llT |. corrected by 
a factor of 1.007 for QED effects, we obtain \V c b\ = 
(39.8 ± 1.8 ± 1.3 ± 0.9 F f) x 10~ 3 , where the third er- 
ror is due to the theoretical uncertainty in <?(1). As 
an alternative, we use a quenched lattice calculation 
based on the Step Scaling Method (SSM) [ijj], and obtain 
|V C J = (41.6 ± 1.8 ± 1.4 ± 0.7 FF ) x 10~ 3 . The authors 
of [12| report the lattice determination of Q(w) for finite 
momentum transfer. Although quenched, this new calcu- 
lation allows the extraction of \V c b\ with relatively small 
model dependence avoiding the large extrapolation to 
w = 1. For example, from a linear interpolation around 
w = 1.2, we obtain IKjI = (41.4±1.3±1.4±1.0 FF )x 10~ 3 . 
We report our measurements of £/(tu)|14;,| for w > 1 in 
Ref. [23|. 

The resulting value of \V c b\ extrapolated to w = 1 is 
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largely independent of previous BABAR results [4[ and 
significantly more precise than previous measurements 
from B — > D£~9( decays. It is also consistent with the 
measurement obtained from B — > D*£~i>£ decays and 
with the inclusive determination of \V c b\ — (41.6 ±0.6) x 
10" 3 [|. 
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TABLE II: Systematic uncertainties in the measurement of <5(1)|L C (,|, p 2 and the branching fraction for B — > D£ vi decays. 
We report the relative error (in %) for <5(l)|C c b and branching fraction, and the absolute error on p 2 . 



Systematic uncertainty on 

| Veil, P 2 and BF 













|C c6 |(%) 


P 2 


BF (%) 


\Vcb\{%) 


P 2 


BF (%) 


|Kb|(%) 


P 2 


BF (%) 


Tracking efficiency 


0.5 


0.008 


0.7 


1.1 


0.003 


1.4 


0.7 


0.004 


1.0 


Neutral reconstruction 


1.0 


0.003 


1.2 


0.8 


0.006 


0.9 


0.9 


0.004 


1.2 


Lepton ID 


1.0 


0.009 


1.0 


0.9 


0.009 


0.8 


0.9 


0.009 


0.9 


Final State Radiation 


0.1 


0.005 


0.2 


0.1 


0.005 


0.2 


0.1 


0.005 


0.2 


Cascade B — > X — > l~ decay background 


0.6 




1.2 


1.0 




2.0 


0.8 




1.5 


B° - B ± cross-feed 


0.2 


0.003 


0.2 


0.2 


0.003 


0.2 


0.2 


0.003 


0.2 


B — ► D*l~i>e form factors 


0.6 


0.008 


0.5 


0.2 


0.003 


0.2 


0.4 


0.006 


0.3 


B — > D**£~i>e form factors 


0.2 


0.007 


0.2 


0.3 


0.006 


0.2 


0.3 


0.007 


0.1 


D branching fractions 


1.0 




2.0 


1.4 




2.7 


1.1 




2.2 


B(B -> D'Tvt) 


1.2 


0.023 


0.6 


1.0 


0.011 


0.9 


1.1 


0.019 


0.6 


B{B — > X£~v t ) 


0.9 




1.9 


0.9 




1.9 


0.8 




1.7 


Stag selection 


1.1 


0.021 


0.6 


1.8 


0.036 


0.8 


1.5 


0.028 


0.8 


B X£~D e fit 


0.7 




1.4 


1.1 




2.2 


0.8 




1.7 


B Dl'vt fit 


1.3 


0.018 


1.1 


1.1 


0.027 


0.6 


1.3 


0.020 


0.8 


B meson lifetime 






0.7 






0.6 






0.6 


Total systematic error 


3.1 


0.04 


4.1 


3.6 


0.05 


5.0 


3.3 


0.04 


4.3 



TABLE III: Fit results for |F c i ) |5(«') extracted at different value of w. In order to reduce any fit model dependence, we fit 
the data interpolating few bins around w. In particular we use 4 bins between w = 1.00 and w — 1.24 to extract |F c t|CJ(w) at 
w — 1.03, 1.05 and 1.10, and 4 bins between w — 1.06 and w — 1.30 to extract |F c [,|C/(io = 1.20). The statistical correlation 
between the measurement at w = 1.20 and the others is p CO n = 0.57, instead the systematic uncertainties can be assumed 
correlated at 100%. In the last column we report the results for |V c f,| obtained using the results for the Q(w) computed in G.M. 
de Divitiis et al, Phys. Lett. B655, 45 (2007). 



w 


\v cb \-g{w)-w 


Q{w) 


|F c6 |-10 J 


1.03 


40.9L5.7L1.3 


1.001L0.019 


40.9±5.7±1.3±0.8 


1.05 


40.2L5.0L1.3 


0.987±0.015 


40.7±5.1±1.3±0.6 


1.10 


38.3L3.3L1.3 


0.943±0.011 


40.6±3.5±1.4±0.5 


1.20 


35.3L1.1L1.2 


0.853L0.021 


41.4±1.3±1.4±1.0 



